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Reply 


J. T. Gosling 1 , P. Riley 2 , and R. M. Skoug 1 

1. Los Alamos National Laboratory, Los Alamos, New Mexico 

2. Science Applications International Corporation, San Diego, California 

We strongly disagree with the essence of the Osherovich [this issue] (hereafter Osherovich) 
comment on one of our papers [Gosling, 1999]. The following paragraphs provide the basis of 
our disagreement and elaborate on why we believe that none of the concluding statements in his 
Comment are true. Our most important point is that one can apply the model developed by 
Osherovich and colleagues to real data obtained at a single point in space to determine the 
polytropic index within magnetic clouds if and only if the highly idealized assumptions of that 
model conform to physical reality. There is good reason to believe that those assumptions do not 
provide an accurate physical description of real magnetic clouds in the spherically expanding solar 
wind. 

The polytropic equation that relates pressure, P, and mass density, p, in a parcel of solar wind 
plasma as a function of time and distance has the form 

Pp - "^ = constant, (I) 

where y is the polytropic index. When one also assumes an ideal gas, this takes the form 

T = Sn Y ~ 1 , (2) 


where T is the temperature, n is the number density, and S is the entropy. Rigorously, to 
determine y it is necessary to sample temperature and density within the same plasma parcel at 
different radial distances from the Sun. In practice, the requirement of sampling the same plasma 
parcel can be relaxed if different plasma parcels sampled at different heliocentric distances start 
with essentially the same entropy and the entropy remains constant with distance as, for example, 
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appears to be the case for protons in the high-speed wind at high heliographic latitudes near solar 
activity minimum [Feldman etal., 1998]. 


Osherovich et al. [1993a] model magnetic clouds in the solar wind as cylindrical magnetic flux 
ropes with axially symmetric plasma and magnetic field properties that are invariant along extended 
cylindrical surfaces. In such a model the function = rA^ is a constant on any cylindrical surface 



azimuthal component of the vector magnetic potential (in cylindrical coordinates). As noted in his 
comment, Osherovich and colleagues also assume a polytropic relationship between gas pressure 
and mass density, an ideal gas equation of state, and uniform cylindrical expansion about the flux 
rope axis as a function of time. If such a model describes physical reality, then Osherovich is 
correct in stating that T and p on expanding cylindrical Surfaces within the flux rope should be 
related to one another by 


T = F('F)p Y ' 1 , 


(3) 


where F( V F) is a constant for any given expanding cylindrical surface, but may vary from one 
surface to another. This equation simply states that all plasma parcels on an expanding cylindrical 
surface are equivalent in that they evolve in time in an identical fashion, as is required by the 
symmetry assumptions of the model. Since Osherovich and colleagues have assumed an ideal gas 
and a polytropic relationship between gas pressure and mass density, (2) must be satisfied as well, 
so that F( V F) is directly related to entropy by 

F(T / ) = S( v F)/m Y ' 1 , (4) 

where m is the average particle mass. They do not seem to have appreciated this simple 
connection between F^F) and entropy; their introduction of FOP) as something different from 
entropy has simply confused the issue. In reality, (3) reduces to (2) on expanding cylindrical 
surfaces and therefore does not depend explicitly upon a flux rope magnetic field topology; the 
essential assumptions are axial symmetry, invariant plasma properties along extended cylindrical 
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surfaces, an ideal gas, a polytropic relationship between P and p, and uniform cylindrical 
expansion. If a flux rope of this nature were to pass over a spacecraft, each expanding cylindrical 
surface would be crossed twice, so that one would obtain pairs of temperature and density 
measurements of essentially equivalent plasma parcels on each expanding cylindrical surface 
within the flux rope. From these paired measurements one could make a series of two-point 
determinations of both FOP) and y, as Osherovich claims. 


Spacecraft measurements indicate that electron temperature, T e , and density, n e , are often 
negatively correlated as a magnetic cloud passes over a spacecraft. Osherovich and colleagues 
[e.g., Osherovich et al, 1993 a,b, 1995, 1998, 1999 \ Fainbergetd., 1996; Farrugiaet al., 1995, 
1999; Osherovich and Burlaga, 1997] have used such measurements together with their model 
assumptions to infer that the electron polytropic index, y>, is less than 1 .0 within magnetic clouds 
and thatF('P) typically is roughly constant in a given flux rope; in addition, they infer the presence 
of multiple flux ropes in some magnetic clouds. For different flux ropes they find different values 
of FOP) and different values of y e (but always < 1 .0). If these determinations of y e are valid, 
then T P must increase as the clouds expand and S( X P) must be nearly constant throughout a given 
flux rope. 

We wish to emphasize that the technique used by Osherovich and colleagues to determine y e 
within magnetic clouds is valid if and only if their model assumptions accurately describe the 
physical properties of real magnetic clouds. If magnetic clouds are not axially symmetric with 
uniform properties along extended cylindrical surfaces and do not expand in a purely cylindrical 
sense, then one does not obtain a series of 2-point samplings of essentially equivalent plasma 
parcels as a cloud passes over a spacecraft. The onus is on Osherovich and colleagues to prove 
that their model assumptions conform to physical reality and apply to real magnetic clouds in the 
spherically expanding solar wind. In this regard, the observed negative correlation between Tg 
and n e within magnetic clouds is not proof that their assumptions are valid or a demonstration 
either that y e is less than 1.0 in magnetic clouds or thatS is constant throughout a given cloud. 
Negative correlations between temperature and density at a single point in the heliosphere can and 
do arise for other reasons [e. g.. Gosling, 1999]. We think it is highly unlikely that real flux 
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ropes in the solar wind have axially symmetric and uniform properties along extended cylindrical 
surfaces, expand in a purely cylindrical sense, and have constant entropy throughout. 

Osherovich notes that in their model the plasma density on a given expanding cylindrical surface 
must decrease with time as the flux rope expands. He also claims that the magnetic cloud 
observed by Ulysses on June 10, 1993 at 4.64 AU contained two separate flux ropes that passed 
Ulysses on June 10.00 - June 1 1.75 and June 1 1.75 - June 13.0, respectively. Figure 1 shows a 
plot of the proton density observed by Ulysses for these and the surrounding time intervals. 
Closest approach to the axis of the first flux rope or “tube” would have occurred approximately on 
June 10.88 and that of the second on June 12.38. It is immediately apparent that for many of the 
pairs of points equidistant in time on either side of closest approach to the tube axes, which would 
correspond roughly to the same expanding cylindrical surfaces in their model, the densities are not 
lower at the later times, in contrast to the model’s assumption. The same is true if one uses the 
observed flow speed to convert time into a spatial distance. Clearly, these presumed separate flux 
rope tubes do not have uniform and axially symmetric plasma properties. Yet, the June 10, 1993 
Ulysses event is often touted as a prime example supporting the model [e.g., Fainberg et al., 
1996; Osherovich and Burlaga, 1997; Osherovich et al., 1999; Osherovich, this issue]. 

Making the further assumption that flux rope expansions are “self-similar”, Osherovich and 
colleagues have found that expansion does not occur unless y e < 1.0 [e.g., Osherovich et al, 
1993 b\ Osherovich, this issue]. Since most clouds are observed to expand, he thus finds support 
for their model in their presumed experimental determinations of y e within magnetic clouds. But, 
as we have noted, those determinations of y e are themselves strongly dependent on the model 
assumptions. One cannot prove the validity of a model in this way. Further, flux rope expansion 
in the solar wind has been demonstrated in one-fluid MHD simulations where y has explicitly 
been chosen to be 1.67 [Vandas etal., 1996a, b], as well as in 3D analytical MHD calculations for 
values of y ranging from 1.2 to 1.67 [Chen and Garren, 1993; Chen, 1996, 1997], Chen[ 1996], 
who does not assume that flux ropes are axially symmetric cylinders, ascribes this difference in 
results to the difference between a 3D and a 2D calculation and Osherovich et al.’s neglect of the 
radial variation in pressure in the background solar wind in which real magnetic clouds are 
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embedded. 


Our own one-fluid simulations have been concerned primarily with dynamic effects associated 
with coronal mass ejection, CME, expansions, such as the forward-reverse shock pairs associated 
with expansions driven by high internal CME pressures [Gosling etal. , 1994 a,b, 1998; Gosling 
and Riley, 1996; Riley et al., 1997]. Usually in these and similar simulations [e.g., Odstrcil and 
Pizzo, 1999a, b] it is assumed thaty = 1.67. However, as we have previously noted [e.g.. 
Gosling, 1999], other choices of y > 1 .0 have only a minor effect on the overall dynamics of 
these expansions, which is where our primary interest lies in these particular simulations. We 
have never claimed that a choice of y = 1 .67 provides a reasonable estimate of electron temperature 
evolution within CMEs. Indeed, we believe that smaller (but > 1 .0) values probably better 
describe the electron thermal evolution within expanding CMEs in the solar wind. Thus the 
comparison shown in Figure la by Osherovich, which shows predictions for y = 1 .67 and which 
Osherovich claims is what we would predict, is not what we would advocate for electrons within 
CMEs in the solar wind, be they magnetic clouds or otherwise. Moreover, in attempting to make 
the comparison seem as poor as possible Osherovich conveniently ignores the fact that electron 
temperatures within CMEs in the solar wind at 1 AU often are considerably greater than 1 x 10 5 K 
[Gosling et al., 1987], We find no basis for his claim that “...no polytropic index can 
accommodate the evolution of temperature in magnetic clouds, if Formula (4) is used.” We also 
stand by our previous comments [Gosling, 1999] to the effect that there is no observational 
evidence to support the idea that electron temperatures within expanding magnetic clouds, or 
expanding CMEs in general, increase with increasing heliocentric distance at any distance from the 
Sun. Indeed, the evidence is to the contrary. For example, we have recently obtained widely 
separated (in heliocentric distance) 2-point measurements of the same magnetic cloud and found 
that both T« and n~ decreased with increasing heliocentric distance despite the fact that at each 
spacecraft T e and n e were negatively correlated [Skoug et al., 2000b]. 

In order to circumvent some of the difficulties associated with predictions of ever increasing 
values of electron temperature as a magnetic cloud expands, Osherovich now incorporates an 
electron heat flux term into the polytropic equation, where he seems to advocate that the only heat 
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flux that affects electron temperature within a flux rope is through the walls of the flux rope tube, 
transverse to the magnetic field. We fail to understand how electrons can carry a significant flow 
of heat transverse to the magnetic field in the essentially collisionless solar wind or how this new 
version of their model makes any more physical sense than the previous one. 

Osherovich claims that the choice of y > 1 .0 in our simulations leads to a positive correlation 
between T and p at a fixed point in space, which he notes is contrary to electron observations 
within magnetic clouds. However, in many of our simulations we obtain neither a linear nor a 
positive correlation between log T and log p (or n) at a fixed point in space [Riley et al, 2000], 
The relationship between T and n at a fixed point depends on the nature of the initial perturbation at 
the inner boundary of the simulation and on heliocentric distance. We obtain a straight line with 
positive slope in log T vs. log n space at a fixed point in space only when the plasma introduced at 
the inner boundary is isentropic and remains unshocked as it propagates outward from the Sun. 
When both of these conditions are satisfied the slope of the log T vs. log n relationship at a single 
point does provide a direct measure of the value of y used in the simulation. 

In some of our simulations we obtain strong negative correlations between T and n at a single 
point in space [Riley et al., 2000]. The upper panel of Figure 2 displays a random set of points in 
log T vs. log n space used to initialize a disturbance lasting 10 hours at 0.13 AU in a spherically 
symmetric, 1 -dimensional, 1 -fluid simulation in which y= 1.5. The speed at 0.13 AU was held 
constant at 702 km s'* and the random density/temperature disturbance was preceded and 
followed by extended intervals of constant density and temperature, and thus also constant 
entropy. The cross marks the value of density and temperature for this ambient surrounding wind. 
The simulation follows the evolution of the disturbance out into the heliosphere; the lower panel of 
Figure 2 shows the resulting log T vs. log n plot at 3.0 AU where the disturbance has evolved into 
a structure that is nearly in pressure balance. The cross in this panel denotes undisturbed ambient 
solar wind at 3 AU. All values of density and temperature are lower at 3.0 AU than at 0.1 3 AU 
owing to the spherical expansion of the solar wind and our choice of y > 1.0. Points in the 
perturbed ambient wind (small dots) lie roughly along a straight line with positive slope (0.5), 
indicating that the perturbation to the surrounding, initially isentropic, ambient wind produced by 
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the initial disturbance was devoid of strong shocks. On the other hand, points within the original 
disturbance pulse (large dots) lie roughly along a line of negative slope (-0.908) even though y for 
the entire simulation was 1.5. The negative slope is a consequence of non-isentropic structure in 
the initial disturbance and evolution toward pressure balance. Thus, negative correlations between 
density and temperature at a single point in space within a CME or magnetic cloud may simply 
reflect the presence of (non-isentropic) structure in the initial disturbance close to the Sun and the 
tendency for such structure to achieve pressure balance as it evolves outward [Gosling, 1999; 
Riley etal., 2000 -,Skoug etal., 2000a]. And, of course, one cannot directly infer the value of y 
from the log T vs. log n plot within the simulated CME in this case. 

Finally, Osherovich and colleagues give the impression that negative correlations between T e and 

n are not observed in the solar wind except within magnetic clouds. However, recent work 

c 

[Skoug etal, 2000a] demonstrates that negative correlations between T e and n e are common in 
much of the solar wind, including within CMEs in the solar wind that are not magnetic clouds. 
This is consistent with previous studies of solar wind thermal core electrons [Phillips and Gosling, 
1990; Hammond et al, 1996]. That is, the negative T e -n e correlation in the solar wind does not 
depend on the presence of a magnetic structure with a flux rope topology. It has also been 
demonstrated that at 1 AU halo electrons do not contribute more to the total pressure in magnetic 
clouds than they do in other types of solar wind of comparable density [Skoug et al., 2000a], 
contrary to statements by Osherovich. 

To summarize, Osherovich’s claim that essentially two-point, rather than single-point, 
measurements of T e and n e are made as a magnetic cloud passes over a spacecraft is valid if and 
only if the the highly idealized symmetry and uniformity assumptions of their model correspond to 
reality. If those assumptions are not physically correct then such measurements reveal nothing 
about the value of y e within magnetic clouds. It is up to Osherovich and his colleagues to show 
that their assumptions conform to physical reality for real magnetic clouds in the spherically 
expanding solar wind, but we have provided solid evidence that they do not. We have also 
emphasized: (1) Their F('F) is simply related to entropy. (2) A value of y e less than 1.0 is not 
required for flux rope expansion in the solar wind. (3) T and p are often not positively correlated 
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in simulations where y > 1 .0. (4) Strong negative correlations between T e and n e at a single point 
in space can arise from non-isentropic structure within magnetic clouds and the solar wind in 
general even when y e > 1 .0. (5) Negative correlations between T e and n e are observed in much of 
the non-magnetic cloud solar wind, making it unlikely that those correlations are related to the 
special magnetic topology of magnetic clouds or are a consequence of y e < 1.0. (6) Available 
evidence indicates that electron temperatures within magnetic clouds, and within CMEs in general, 
decline as those disturbances expand out into the heliosphere. (7) If a polytrope equation 
adequately describes electron evolution, the polytropic index is greater than 1.0 in magnetic 
clouds, in non-cloud CMEs, and in the solar wind in general. 
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Figure Captions 


Figure 1 . Solar wind density measured by Ulysses at 4.64 AU and S32.5° during passage of a 
disturbance driven by an overexpanding CME that was also a magnetic cloud. The two separate, 
but adjacent, flux ropes within the CME/cloud identified by Osherovich and colleagues are 
indicated by the solid vertical lines. Dashed vertical lines indicate closest approach to the 
presumed axes of those flux ropes. 

Figure 2 . (a) Density-temperature points used to initialize a 10-hour disturbance in the solar 
wind at 0. 13 AU in a one-dimensional (radial), one-fluid simulation. The cross marks the value of 
density and temperature in the surrounding ambient solar wind at 0. 1 3 AU. (b) The resulting 
paired values of density and temperature at 3.0 AU. The solid curve is the best-fit straight line to 
plasma parcels that originated within the 10-hour disturbance at 0.13 AU. See text for full 
explanation. Adapted from Riley et al. [2000]. 
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